AFRL-RX-WP-JA-2018-0086 


PREDICTING  FATIGUE  DAMAGE  OF  COMPOSITES 
USING  STRENGTH  DEGRADATION  AND 
CUMULATIVE  DAMAGE  MODEL  (POSTPRINT) 


Arafat  I.  Khan  and  Satchi  Venkataraman 
San  Diego  State  University 


Ian  Miller 

N&R  Engineering 


22  SEPTEMBER  2017 
Interim  Report 


Distribution  Statement  A. 

Approved  for  public  release:  distribution  unlimited. 


©  2018  MDPI 


(STINFO  COPY) 


AIR  FORCE  RESEARCH  LABORATORY 
MATERIALS  AND  MANUFACTURING  DIRECTORATE 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7750 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


REPORT  DOCUMENTATION  PAGE  OMBNoXmfws 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 
Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YY)  2.  REPORT  TYPE 

22  September  2017  Interim 

3.  DATES  COVERED  (From  -  To) 

24  May  2016-22  August  2017 

4.  TITLE  AND  SUBTITLE 

PREDICTING  FATIGUE  DAMAGE  OF  COMPOSITES  USING 
STRENGTH  DEGRADATION  AND  CUMULATIVE  DAMAGE 
MODEL  (POSTPRINT) 

5a.  CONTRACT  NUMBER 

FA8650-16-C-5010 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

65502F 

6.  AUTHOR(S) 

1)  Arafat  I.  Khan  and  2)  Ian  Miller  - 

Satchi  Venkataraman  -  N&R  Engineering 

San  Diego  State  University 

5d.  PROJECT  NUMBER 

3005 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

X155 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

1)  San  Diego  State  University  2)  N&R  Engineering 

5500  Campanile  Drive  6659  Pearl  Rd  #  201 

San  Diego,  CA  92182  Cleveland,  OH  44130 

(continued  on  page  2) 

8.  PERFORMING  ORGANIZATION  REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 

Materials  and  Manufacturing  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7750 

Air  Force  Materiel  Command 

United  States  Air  Force 

10.  SPONSORING/MONITORING  AGENCY 
ACRONYM(S) 

AFRL/RXCM 

11.  SPONSORING/MONITORING  AGENCY 

REPORT  NUMBER(S) 

AFRL-RX-WP-JA-20 18-0086 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Distribution  Statement  A.  Approved  for  public  release:  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 

PA  Case  Number:  8 8ABW-20 17-4652;  Clearance  Date:  22  Sep  2017.  This  document  contains  color.  .  Journal  article  published 
in  merican  Journal  of  Composites  Science,  Vol.  2,  No.  9,  14  Feb  2018.  ©  2018  MDPI.  The  U.S.  Government  is  joint  author  of 
the  work  and  has  the  right  to  use,  modify,  reproduce,  release,  perform,  display,  or  disclose  the  work. 

The  final  publication  is  available  at  www.mdpi.com/journal/jcs  doi:10.3390/jcs2010009 _ 

14.  ABSTRACT  (Maximum  200  words) 

This  paper  discusses  implementation  of  an  empirical  strength  and  stiffness  degradation  scheme  to  predict  fatigue  failure  and 
residual  stiffness  and  strength  of  laminated  composites.  A  cumulative  damage  accumulation  approach  was  implemented  to 
account  for  the  non-constant  stress  levels  that  occur  due  to  stress  distributions  arising  from  stiffness  degradation  and  local 
failures  during  fatigue  loading.  Degradation  of  strength  and  stiffness  during  fatigue  loading  of  the  composite  was  implemented 
following  the  empirical  model  presented  by  Shokrieh  and  Lessard.  The  fatigue  analysis  was  performed  with  AB  AQUS™  finite 
element  software  employing  a  user-defined  material  subroutine  (UMAT)  developed  for  the  material  response.  Implementation 
results  were  first  verified  for  unidirectional  laminate  test  cases,  validated  by  predicting  stress  versus  life  (S-N)  curves  for 
several  laminate  coupons  test  simulations  and  residual  strengths  of  Open  Hole  Tension  (OHT)  specimens  subjected  to  constant 
amplitude  fatigue  loading. _ 


15.  SUBJECT  TERMS 

Empirical  strength  and  stiffness  degradation;  fatigue  failure;  residual  stiffness;  laminated  composite;  user-defined  material 
subroutine  (UMAT);  unidirectional  laminate  test  case;  stress  versus  life  (S-N)  curve;  Open  Hole  Tension  (OHT);  amplitude 
fatigue  loading 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT: 

SAR 

18.  NUMBER 

OF  PAGES 

23 

19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 
Jeff  Simmons 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

(937)  255-1346 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


composites  science 


MDPI 


Article 

Predicting  Fatigue  Damage  of  Composites  Using 
Strength  Degradation  and  Cumulative 
Damage  Model 

Arafat  I.  Khan  lf*,  Satchi  Venkataraman  1  B  and  Ian  Miller  2 

1  Department  of  Aerospace  Engineering,  San  Diego  State  University,  San  Diego,  CA  92182,  USA; 
satchi@mail.sdsu.edu 

2  N  &  R  Engineering,  6659  Pearl  Rd  #  201,  Cleveland,  OH  44130,  USA;  imiller@nrengineering.com 
*  Correspondence:  arafat.khan@sdsu.edu;  Tel.:  +1-619-594-4777 

Received:  30  December  2017;  Accepted:  8  February  2018;  Published:  14  February  2018 

Abstract:  This  paper  discusses  the  prediction  of  fatigue  response  of  composites  using  an  empirical 
strength  and  stiffness  degradation  scheme  coupled  to  a  cumulative  damage  accumulation  approach. 
The  cumulative  damage  accumulation  approach  is  needed  to  account  for  the  non-constant  stress 
levels  that  arise  due  to  stress  distributions  from  stiffness  degradation  during  the  fatigue  loading. 
Degradation  of  strength  and  stiffness  during  fatigue  loading  of  the  composite  was  implemented 
by  following  the  empirical  model  presented  by  Shokrieh  and  Lessard  with  some  modification  and 
correction  to  the  non-dimensional  load  parameter  definition.  The  fatigue  analysis  was  performed 
using  ABAQUS™  finite  element  software  using  a  user-defined  material  subroutine  UMAT  developed 
for  the  material  response.  Implementation  results  were  first  verified  for  unidirectional  laminate  test 
cases  and  validated  by  predicting  stress  versus  life  (S-N)  curves  for  several  laminate  coupons  test 
simulations  and  residual  strengths  of  Open  Hole  Tension  (OHT)  specimens  subjected  to  constant 
amplitude  fatigue  loading. 

Keywords:  CFRP  composite;  fatigue;  strength;  stiffness;  UMAT 


I.  Introduction 

The  Composite  materials  used  in  primary  aircraft  structures  can  experience  fatigue  damage  and 
failure  due  to  the  repeated  loads  they  experience.  Tools  and  models  that  stimulate  the  response  of 
composites  under  cyclic  loads  are  necessary  in  order  to  design  structures  that  can  withstand  such 
damage.  Effective  tools  can  also  help  to  estimate  remaining  lifetimes  and  residual  strength  and  stiffness 
of  composites  as  they  are  subjected  to  repeated  or  cyclic  loading. 

The  response  of  composite  structures  under  fatigue  loading  investigated  by  a  number  of  researchers 
that  has  led  to  development  of  a  number  of  fatigue  prediction  models.  These  models  can  be  classified  in  the 
following  categories  [1]:  fatigue  life  concepts  [2,3],  phenomenological  models  with  stiffness  [4,5],  strength 
degradation  models  [6,7],  continuum  damage  mechanics  (CDM)  based  models  [8-10]  and  micromechanics 
models  [11,12]  and  uncertainty  and  Bayesian  based  probabilistic  models  [13-15]. 

A  commonly  used  approach  in  fatigue  life  predictions  is  to  use  stress  versus  life,  known  as  S-N 
curves.  These  fatigue  life  models  [2,3]  use  test  data  from  constant  amplitude  tests  at  varying  stress 
levels  and  fit  empirical  S-N  curves  that  relate  the  number  of  cycles  to  failure  to  the  applied  stress 
level.  The  constant  amplitude  cyclic  loads  are  characterized  by  the  mean  stress  level  (crm)  and  the 
amplitude  (cra)  of  the  stress  variations  about  the  mean.  This  is  alternatively  expressed  in  terms  of 
the  maximum  stress  ( cmax  =  crm  +  cra)  and  the  stress  ratio  or  R-ratio  (R  =  amin / crmax) •  Goodman  type 
diagrams  [3]  provide  models  that  can  fit  the  cycles  to  failure  as  a  function  of  mean  stress  and  varying 
stress  amplitude.  This  requires  extensive  experimental  work,  as  measuring  the  cycles  to  failure  for  all 
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possible  mean  and  alternating  stress  values  states  that  do  not  lead  to  static  stress  failure  is  quite  large. 
The  laminates  chosen  for  such  tests  should  be  able  to  capture  the  various  modes  of  intra-laminar  (fiber 
and  matrix)  failure  within  a  ply.  A  limitation  of  this  approach  is  that  at  low  stress  levels,  it  is  expensive 
to  perform  the  cyclic  tests  until  the  specimen  fails.  As  a  result,  the  tests  are  truncated  at  number  of 
cycles  that  represents  an  upper  limit  of  what  an  engineering  component  would  experience  in  service. 

Fatigue  response  of  composite  materials  can  be  described  by  the  stiffness  and  strength  degradations 
they  exhibit  due  to  accumulation  of  damage.  These  are  characterized  by  the  residual  stiffness  and  strength 
measurements  of  components  tested  at  different  stress  levels  and  to  varying  fractions  of  their  life  (previously 
determined  for  life  models).  Phenomenological  models  based  on  empirical  fits  to  the  test  data  provide 
evolution  laws  that  describe  a  gradual  reduction  in  the  laminate  stiffness  and  strength  at  macroscopic 
scale  due  to  the  damage  evolution  at  microscopic  constituent  length  scales.  A  number  of  models  have 
been  proposed  for  stiffness  degradation  during  fatigue  [4,5]  and  strength  degradation  [6,7].  In  strength 
degradation  models,  the  ply  fails  when  the  strength  value  reaches  the  ply  stress  level. 

Fatigue  behavior  in  composites  has  also  been  modelled  using  continuum  damage 
mechanics-based  approaches  by  using  fatigue  cycles  to  account  for  damage  growth  [8-10].  These  types 
of  approaches  suggest  that  degradation  of  the  composite  material  is  in  direct  relation  with  specific 
damage  such  as  matrix  cracking  or  delamination.  The  continuum  damage  mechanics  (CDM)  models 
correlate  one  or  more  material  properties  with  defined  damage  variables  which  are  determined 
quantitatively  from  the  damage  progression  or  failure  mechanism. 

One  of  the  new  category  of  approaches  to  model  fatigue  in  composites  are  models  based  on  uncertainty 
and  Bayesian  based  probabilistic  framework  [13-15].  Most  of  the  probabilistic  models  are  based  on 
quantifying  model  uncertainties  for  different  fatigue  models.  A  progressive  damage  in  composite  material 
using  a  full  Bayesian  prediction  and  updating  framework  is  presented  by  Chiachio  et  al.  [13].  They  used 
a  cumulative  damage  model  based  on  the  theory  of  Markov  chains  to  infer  a  complete  damage  progression. 
Chiachio  et  al.  [14]  presented  one  such  model  based  on  the  uncertainty  for  empirical  fatigue  model. 
The  authors  did  not  account  for  any  interlaminar  failure  for  the  Bayesian  approach.  In  the  work  of 
Shankaraman  et  al.  [15],  several  different  models  such  as  finite  element  (FE),  surrogate  model  analysis, 
crack  growth  model  and  the  effect  of  the  types  of  uncertainties  such  as  physical  variability,  data  uncertainty 
and  errors  were  incorporated  in  both  the  inverse  problem  and  the  forward  problem.  Diao  et  al.  [16] 
incorporated  a  statistical  approach  for  reinforced  composite  laminates  subjected  to  multiaxial  loading  for 
the  static  and  fatigue  properties.  They  obtained  the  statistical  distribution  function  of  the  material  strength 
parameters  by  fitting  experimental  data  to  a  two-parameter  Weibull  distribution. 

After  a  review  of  material  models  proposed  in  the  literature,  we  selected  an  empirical  model  proposed 
by  Shokrieh  and  Lessard  [6,17-19].  Our  model  tests  six  basic  unidirectional  laminates  under  fatigue  load 
to  produce  fatigue  data  for  AS4  3501-6  graphite  epoxy  material.  A  material  user  subroutine  (UMAT)  for 
modeling  the  material  degradation  during  fatigue  loading  was  developed  for  use  within  the  ABAQUS™ 
finite  element  software.  A  cumulative  damage  approach  was  used  to  account  for  the  change  in  stress 
amplitude  resulting  from  the  stress  state  redistribution  after  failure.  The  material  degradation  model 
is  based  on  life  prediction  models  that  use  empirically  fitted  constant  life  curves  [20]  to  determine  the 
number  of  cycles  to  failure  and  a  degradation  model  where  the  reductions  in  stiffness  and  strength  occur 
as  a  function  of  the  number  of  cycles  of  loading  experienced  as  a  fraction  of  the  total  life  time. 

2.  Constant  Life  Diagrams 

Constant  life  diagrams  (CLD)  are  obtained  from  empirical  fits  to  fatigue  test  data  that  relate 
cycles  to  failure  to  the  mean  and  alternating  stress  values  of  the  constant  amplitude  fatigue  loading. 
The  constant  life  curves  are  simply  the  contour  curves  corresponding  to  a  given  number  of  cycles  of 
failure,  plotted  in  the  space  of  mean  and  alternating  stress.  Many  different  models  and  methods  have 
been  proposed  to  parametrically  describe  these  curves  [3].  Figure  la  shows  a  typical  constant  life 
diagram  (CLD).  The  main  axes  on  a  CLD  graph  [20]  are  the  mean  cyclic  stress  (cm)  and  the  cyclic  stress 
amplitude  ( cra ).  In  addition,  radial  lines  of  constant  R- ratio  are  shown  on  the  graph  to  easily  recognize 
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how  the  stress  ratio  affects  the  shape  of  these  curves.  Goodman  diagrams  use  linear  fits  for  constant 
life  curves.  However,  for  composite  materials  researchers  find  that  constant  life  curves  are  nonlinear 
and  non-symmetric  with  respect  to  the  mean  stress  axis.  These  curves  allow  interpolation  of  test  data 
to  estimate  lifetimes  for  given  stress  loadings.  A  number  of  analytical  models  for  parametrizing  the 
curves  have  been  published  in  literature  [3,20-22].  Adam  et  al.  [20]  presented  an  analytical  model  that 
not  only  parametrizes  the  CLD  but  also  allows  unifying  the  different  curves  into  a  single  two-parameter 
fatigue  curve  which  in  this  work  is  referred  to  as  the  normalized  fatigue  life  model  (Figure  lb). 


Figure  1.  (a)  Schematic  of  a  typical  constant  life  diagram;  (b)  Normalized  fatigue  life  model  derived 
from  constant  life  diagram. 


A  two-parameter  model  of  a  constant  life  diagram  was  fitted  to  the  non-dimensional  mean  stress 
q  and  stress  amplitude  a  values  and  the  ratio  c  of  the  static  strengths  Srej\  and  Sref2  of  the  orthotropic 
ply  material  as  shown  in  Equation  (1). 


cra  Sref2  crm 

Q - 'C=C  '  V=Q - 

^ref  1  ^ref  1  ^ref  1 


(i) 


For  the  composite  lamina  in  the  fiber  directions,  the  Srej\  and  Sref2  terms  in  Equation  (1)  take  on 
values  of  the  fiber  tension  and  fiber  compression  strengths  Xj  and  Xq,  respectively.  For  the  transverse 
to  fiber  direction,  the  Srej\  and  Srej 2  terms  in  Equation  (1)  take  on  values  of  the  transverse  to  fiber 
compression  and  transverse  to  fiber  tension  strengths  Yq  and  Yj,  respectively.  Normalization  allows 
the  definition  of  a  unified  load  parameter  u  to  fit  constant  life  diagrams.  The  parametric  model 
proposed  by  Adam  et  al.  [16]  is  given  by: 


a=f(l-q)u(q  +  c)°  (2) 

where/,  u  and  v  are  model  fit  constants.  The  exponents  u  and  v  determine  the  shapes  of  the  left  and 
right  wings  of  a  bell-shaped  CLD  curve.  However,  Gathercole  et  al.  [22]  showed  that  the  degree  of 
curve-shape  asymmetry  was  not  very  great  for  carbon-epoxy  composites;  therefore,  the  assumption 
that  u  and  v  are  equal  is  valid.  Thus,  u  is  given  by: 

U  = _ /w(g//) _  (3) 

ln[(l-q)(q  +  c)\  W 

In  Equation  (3),/  is  a  curve  fitting  parameter  with  a  value  of  1.06. 

The  unified  load  parameter  u  is  used  to  relate  the  loading  to  the  number  of  cycles  to  failure  in  the 
empirical  model.  Shokrieh  used  a  linear  fit  model  to  the  logarithm  of  the  number  for  cycles  to  failure 
log(A/)  as  a  function  of  the  unified  load  parameter  u, 

u  =  v  =  A  +  B  log  Nf  (4) 
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where  A  and  B  are  the  curve  fitting  constants. 

The  number  of  cycles  to  failure  (Ny)  is  computed  for  the  state  of  axial  fiber  direction  loading  using 
stress  (Xu  and  for  the  transverse  to  fiber  direction  loading  using  stress  a 22  by  choosing  the  appropriate 
normalization  stress  in  the  definition  of  the  load  parameter  u,  as  described  above.  Shokrieh  and 
Lessard  also  showed  that  the  unified  load  parameter  for  the  pure  shear  state  of  stress  {<J\2,  (t13/  or  a 23) 
requires  modification.  There  is  no  equivalent  loading  direction  strength  tension  or  compression  when 
loaded  in  shear.  The  unified  load  parameter  for  shear  loading  as  proposed  by  Shokrieh  [6]  is 


u  =  log10 


'  Hz/f)  ' 

Un[(l-q)(l  +  q)]_ 


(5) 


In  the  above  expression  z  =  Cal<Js}iear  and  *7  =  Cshm-r/S,  where  ca  is  the  stress  amplitude,  crshear  is 
the  applied  shear  stress  and  S  is  the  static  shear  strength. 

The  unified  parameter  u,  in  Equation  (3)  is  also  modified  slightly  for  loading  in  the  transverse  to 
fiber  direction  (matrix  failure  dominated  state  of  stress),  a 22 •  Our  implementation  of  the  model  found 
that  the  estimation  of  the  number  of  cycles  to  the  failure  using  the  general  form  of  Equation  (3)  used  by 
Shokrieh  [6]  does  not  provides  a  realistic  prediction  for  the  loading  transverse  to  the  fiber.  The  model 
is  required  to  be  normalized  by  the  compression  strength,  where  the  general  form  of  Equation  (3) 
normalized  using  the  tensile  strength  in  the  matrix  direction.  Implementations  of  Shokrieh's  model  by 
other  authors  found  in  the  literature  does  not  address  this  normalization  approach  for  the  unified  load 
parameter  in  the  matrix  direction  [23-26].  The  unified  parameter  in  the  state  of  stress  a 22  is  normalized 
by  the  strength  parameter  Yq  and  the  unified  parameter,  u  can  be  expressed  as: 


_ ln(a/f) _ 

ln[((crm/Yc)  +  1)((YT/YC)  -  Om/Yc))] 


(6) 


The  experimental  value  of  the  transverse  compressive  strength  Yq  is  significantly  higher  than  the 
tensile  strength  Yj,  where  the  longitudinal  tensile  strength  Xj  is  generally  larger  than  the  compressive 
strength  Xq-  Therefore,  for  the  a 22  state  of  stress,  the  normalization  was  performed  by  the  compressive 
strength  of  Yq,  as  shown  in  Figure  2b. 


Figure  2.  Unified  parameter  determination  for  (a)  fiber  direction;  (b)  matrix  direction  loading. 


Table  1  shows  a  comparison  of  the  computed  values  of  the  unified  load  parameters,  u  and  u  using 
the  model  presented  by  Shokrieh  [6]  given  Equation  (5)  and  the  modification  that  was  made  in  this 
paper  which  is  given  by  Equation  (6).  The  accurate  prediction  of  the  unified  parameter,  u  modification 
is  significant  to  predict  fatigue  life  along  the  matrix  direction  as  failure  generally  initiates  in  the  matrix. 
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Table  1.  A  comparison  to  the  values  of  the  unified  load  parameters. 


Stress  Level 

u  (Shokrieh's  Model) 

Nf  (Shokrieh's  Model) 

u  (Modified) 

Nf  (Modified) 

90%  of  Yt 

-1.2247 

6.8624  x  1CT24 

1.20938 

155.394 

80%  of  Yt 

-1.2344 

5.6997  x  l(r24 

1.33436 

3114.24 

70%  of  Yt 

-1.2671 

2.4325  x  1CT24 

1.47316 

86931.2 

60%  of  Yj 

-1.3282 

5.7220  x  10-25 

1.62986 

3.7281  x  106 

50%  of  Yr 

-1.4151 

7.1256  x  10~26 

1.81064 

2.8482  x  108 

3.  Shokrieh's  Material  Degradation 

Shokrieh  and  Lessard  [17,18]  proposed  a  power  law  model  (Figure  3)  to  represent  the  relationship 
between  the  number  of  cycles  and  the  residual  strength  or  stiffness  of  a  laminate.  This  model  is 
capable  of  capturing  both  the  "sudden  death"  [27]  and  "wear  out"  [28]  modes  for  failure  in  laminates. 
The  "sudden  death"  mode  is  observed  in  fatigue  tests  at  high  stress  levels  whereas  the  "wear  out" 
mode  is  observed  in  fatigue  tests  at  low  stress  levels. 


Figure  3.  The  sudden  death  and  wear  out  model  for  the  power  law  fit  used  in  Shokrieh's  model. 


The  power  law  fit  for  the  degradation  of  the  strength  and  the  stiffness  follows  the  form: 


S*  +  N?  =  1  and  EA  +  =  1 

The  number  of  cycles  to  failure  is  normalized  as  follows: 

N  —  ^(”)-%(°^)  ;  0  <  jy  <  ! 
log(Nf)  -  log( 0.25) 

The  residual  strength  and  stiffness  are  normalized  using: 


(7) 


(8) 


Strength  :  S 


(Jr  ~  <J 
Jo-C' 


;  0  <  S  <  1;  Stiffness  :  E  = 


Er-Ef 


E0  -  Ef 


;  0  <  E  <  1 


(9) 


Thus,  the  residual  strength  and  stiffness  as  a  function  of  the  number  of  cycles,  stress  level  and 
stress  ratio  are  given  by: 


C7>(K,C7,ft) 


1  - 


Eog  n  —  Log  0.25 
Log  Nf  —  Log  0.25 


(<7o  -  c)  +  Jo 


(10) 
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E(k,  a,n) 


1 


/  Log  n  -  Log  0.25 
^  Log  Nf-  Log  0.25 


(ii) 


where  ar  is  the  residual  strength,  ctq  is  static  strength,  a  is  the  stress  level,  n  is  the  number  of  applied 
cycles,  Nf  is  the  number  of  cycles  to  failure  at  the  current  stress  level  (fatigue  life),  oc  and  /3  are  the 
experimental  curve  fitting  parameters,  E  is  the  residual  stiffness,  a  is  the  stress  level,  n  is  the  number 
of  applied  cycles,  E$  is  the  initial  stiffness,  £f  is  the  strain  at  failure  at  applied  stress  a  and  7  and  A  are 
the  experimental  curve  fitting  parameters  [17,18].  The  curve  fitting  parameters  oc,  7  and  A  are  given 
in  Table  2.  These  fitting  parameters  provided  by  Shokrieh  and  Lessard  [6,17-19]  contain  typographical 
errors  where  the  oc  and  parameters  are  switched. 


Table  2.  Fatigue  material  constants  in  Shokrieh  model  for  AS4  3501-6. 


Test 

A 

7 

ef 

oc 

0 

A 

B 

Longitudinal  tension 

14.57 

0.3024 

0.0136 

0.473 

10.03 

1.3689 

0.1097 

Longitudinal  compression 

14.57 

0.3024 

0.0136 

0.025 

49.06 

1.3689 

0.1097 

Transverse  tension 

14.57 

0.1155 

0.0068 

0.1255 

9.628 

0.999 

0.096 

Transverse  compression 

14.57 

0.1155 

0.0068 

0.0011 

67.36 

0.999 

0.096 

In-plane  shear 

0.70 

11 

0.101 

9.11 

0.16 

0.0999 

0.186 

Out-of-plane  shear 

0.70 

11 

0.101 

12 

0.20 

0.0999 

0.111 

The  material  constants  determined  for  the  residual  strength  and  the  stiffness  model  reported  by 
Shokrieh  and  Lessard  contains  typographical  error  in  all  the  reported  publications  by  the  authors. 
Figure  4  shows  the  use  of  the  parameters  oc  and  /3  for  the  strength  degradation  model  of  the  uniaxial 
0°  test  specimen  as  given  by  Equation  (10).  It  is  evident  that  the  values  of  the  parameters  are 
switched  [17,18]  and  causes  significant  discrepancies  in  residual  strength  prediction. 


a> 

1 0.8 
<n 


■O 

’<75  0.6 

Q) 

K 

■o 

(V 

N 

15  0.4 


V 

\ 

\ 


\ 


\ 


-  a  =  0.473  and/?  =10.03 
.  a  =  10.03  and  p  =  0.473  (a,  p  reported  by  Shokrieh  and  Lessard) 
-Model  plotted  by  Shokrieh  and  Lessard 


\ 


0  0.2  0.4  0.6  0.8  1 

Normalized  Number  of  Cycles 

Figure  4.  The  normalized  residual  strength  model  for  the  0°  specimen  tested  by  Shokrieh  and  Lessard 
showing  the  effect  of  the  parameters  oc  and  j 6  obtained  from  the  test  data. 


Cumulative  Damage  Approach 

Number  of  authors  have  implemented  Shokrieh  and  Lessards's  empirical  model  to  predict  fatigue 
behavior  of  composites  [23-26].  The  authors  [23-26]  did  not  account  for  any  damage  accumulation 
when  the  state  of  stress  redistributes  due  to  failure  in  a  material  point.  The  above  models  were 
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developed  using  constant  amplitude  tests.  To  apply  them  to  structural  fatigue  simulations  under 
varying  loads,  a  cumulative  damage  approach  [29]  is  needed.  Variations  in  stress  amplitudes  for 
applied  stress  may  be  due  either  to  the  applied  load  being  time  varying  or  to  stress  distributions  arising 
from  stiffness  changes  caused  by  stiffness  degradation  or  local  material  failure.  The  cumulative  damage 
approach  introduces  a  damage  variable  to  track  damage  accumulations  and  maintains  continuity  of 
the  damage  variable  when  stress  changes  occur.  In  this  way  damage  accumulation  at  a  new  stress 
level  starts  with  the  initial  damage  level  being  what  was  sustained  already  due  to  previous  history  at 
different  stress  levels.  The  cumulative  damage  not  only  account  for  the  damage  due  to  the  stress  state 
redistribution  but  also  accounts  for  the  existing  damage. 

If  a  material  that  degrades  due  to  cyclic  loading  is  subjected  to  uniaxial  fatigue  load  with  maximum 
stress  a,  failure  occurs  when  the  residual  strength  <jr,  reaches  the  applied  stress  level  a.  Using  the 
strength  degradation  model  (Equation  (10)),  the  ratio  of  the  stress  margin  due  to  the  degradation  of 
strength  (ar  —  a)  after  n  cycles,  to  the  maximum  stress  margin  at  start  of  the  fatigue  loading  (cq  —  a)  is 
expressed  as: 


f 

1- 

Zft(4ft) 

/  \ 

InUNf) 

(12) 


where  crr,  is  the  residual  strength,  (ft)  is  the  number  of  cycles  at  applied  load  and  (Nj)  is  the  number  of 
cycles  to  failure.  The  above  ratio  varies  from  unity  at  start  to  zero  at  failure.  In  this  work,  we  adopt 
a  definition  of  damage  based  on  the  reduction  in  material  strength.  The  damage  variable  D  therefore 
takes  the  form: 


i 


(13) 


Figure  5  provides  a  schematic  of  the  cumulative  damage  approach.  The  curves  represent  the 
residual  strength  degradation  as  a  function  of  the  number  of  cycles  for  different  constant  amplitude 
fatigue  loading  stress  levels  (cq  <  a2  <  Connecting  the  end  points  of  each  curve  also  gives  the 
S-N  curve.  In  the  cumulative  damage  approach,  we  maintain  equivalency  of  damage  between  fatigue 
cycles  when  there  is  a  stress  level  change  by  equating  the  residual  strength.  Let  us  assume  D\e  is  the 
damage  accumulated  at  the  end  of  the  elapsed  loading  block  of  n\  fatigue  cycles  at  stress  level  <J\. 
If  the  stress  level  changes  to  <j2  after  the  first  block  of  n\  fatigue  cycles  at  stress  level  cq,  the  damage 
accumulation  starts  of  the  new  block  D2s  must  be  equal  to  D\e.  To  use  a  degradation  model  for  constant 
amplitude  loading,  with  the  new  stress  level  cr2  we  calculate  first  the  equivalent  number  of  cycles  ns2 
of  loading  at  stress  level  a2  that  would  have  created  the  damage  level  D\e.  Using  the  definition  of 
damage  in  Equation  (13),  the  equivalent  cycles  ft|  is  calculated  as  shown  below  in  Equation  (14). 


In  (4Ny2 ) 

4ft|  =  (4ft1)ln(4N/l) 


(14) 


The  remaining  life  at  stress  level  cr2  is  nr2  —  Nj2  —  ft|.  If  this  is  greater  than  zero  then  the  fatigue 
loading  can  continue  and  the  material  degrades  following  the  degradation  law  corresponding  to  stress 
level  cr2.  The  damage  accumulation  in  the  new  block  n2  cycles  at  stress  level  a2  is  calculated  based  on 
total  cycle  count  equal  to  (ft|  +  n2).  For  any  fcth  new  block  of  fatigue  load,  we  calculate  the  equivalent 
elapsed  cycles  as  nsk  and  the  remaining  life  nrp  using  Equation  (15). 

H^fk) 

4 n\  -  (4nfc_i)  (n(4N/*-l)  and  nrk  -  Nfk  -  nk  (15) 
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4.  Failure  Criteria  and  Instant  Material  Degradation 


In  our  work,  we  use  Hashin's  [30]  failure  criteria  to  determine  the  fiber  and  matrix  failures  in 
a  unidirectional  lamina.  Equations  (16)-(19)  summarize  equation  for  the  failure  envelopes  for  fiber 
and  matrix  failure  from  Hashin's  criteria. 

Fiber  failure  in  tension 


faU\2  ^l22  +  ^13  , 

UJ  S?2 


(16) 


Fiber  failure  in  compression 


(17) 


Matrix  failure  in  tension  ( C22  +  033  >  0) 


(^22  +  ^33)  ,  ^23  —  ^22^33  ^12  +  a[ 


Y2 

1 T 


+ 


s2 

D23 


+ 


13 


S2 

d12 


=  1 


(18) 


Matrix  failure  in  compression  (C22  +  033  <  0) 


*lV-1 

2  S23J 


C22  +  ^33 


(C722  +  C733  )  2  (^22  +  ^33)2  .  ^23  “  ^22^33  .  ^12  +  °\ 


4S2 

*D23 


y2 

iT 


+ 


s2 

^23 


13 


S2 

d12 


=  1  (19) 


In  Equations  (16)-(19),  Xj  and  Xq  are  the  longitudinal  tensile  and  compressive  strengths,  Yj  and 
yc  are  the  transverse  tensile  and  compressive  strengths,  S\2  is  the  in-plane  shear  strength  and  S23  is 
the  out  of  plane  shear  strength.  An  instantaneous  ply  stiffness  degradation  scheme  is  used  for  the 
progressive  failure  when  ply  failure  is  detected.  Camanho  [31]  showed  that  a  small  but  finite  residual 
stiffness  after  failure  is  needed  to  accurately  predict  the  evolution  of  progressive  failure.  Table  3  shows 
the  residual  stiffness  of  the  ply  following  failure  in  each  mode. 


Table  3.  Material  properties  to  be  degraded  based  on  the  mode  of  failure. 


Failure  Mode 

Ei 

E2 

E3 

G12 

Gl3 

C23 

Fiber  tension 

0.07 

- 

- 

Fiber  compression 

0.14 

- 

- 

Matrix  tension 

- 

0.2 

0.2 

0.2 

0.2 

0.2 

Matrix  compression 

- 

0.4 

0.4 

0.4 

0.4 

0.4 
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5.  Finite  Element  Implementation 

Finite  element  analysis  software  ABAQUS  [32]  was  used  to  perform  the  fatigue  analysis  of 
composite  laminates  and  laminates  structures.  Fatigue  material  degradation  was  implemented 
as  a  user-defined  material  subroutine  (UMAT)  for  use  in  ABAQUS  non-linear  implicit  solution. 
The  laminates  were  modelled  as  a  three-dimensional  solid  with  ply  by  ply  representation.  Each  ply 
was  modeled  using  the  eight-node  bilinear  reduced  integration  C3D8R  solid  elements,  with  enhanced 
hourglass  control.  Figure  6  shows  a  flow  chart  of  the  material  degradation  model  implementation. 
The  Jacobian  update  implemented  in  the  UMAT  subroutine  was  consistent  for  use  with  geometric 
nonlinear  analyses.  Failure  indices,  residual  strengths  and  damage  parameters  at  each  integration 
point  are  stored  as  Solution  Dependent  Variables  (SDVs)  in  ABAQUS. 

The  fatigue  loading  is  simulated  by  two  load  steps:  a  static  load  step  under  load  control  up  to  the 
maximum  loading  during  each  cycle,  followed  by  a  constant  load  over  a  specified  time.  The  degradation 
with  cycles  is  performed  with  a  forward  time  marching  scheme  in  with  an  increment  of  1000  cycles  for 
each  time  increment  in  the  second  step.  To  determine  residual  stiffness,  additional  load  steps  are  added  to 
unload  the  structure  and  then  load  the  structure  under  displacement  control  until  failure. 


Figure  6.  Flow  chart  of  the  FEM  implementation  using  UMAT  subroutine. 

6.  Verification  and  Validation 

The  UMAT  subroutine  implemented  in  ABAQUS  was  first  verified  using  experimental  data  for 
the  six  unidirectional  laminate  cases  presented  by  other  authors  for  stress  ratios  of  R  =  0.1  and  10. 
This  was  followed  by  simulations  of  laminate  coupons  tested  under  homogenous  stress  states.  Figure  7 
shows  the  comparison  of  the  model  predicted  S-N  curve  with  experimental  data  and  simulations 
results  obtained  from  the  literature  for  [30°]i6  laminate  from  Ref.  [6],  [90° ]s  laminate  from  Ref.  [23], 
[0°/90°/90°]s  from  Ref.  [23]  and  [45°/— 45°]2S  from  Ref.  [33].  The  test  data  exhibits  significant  scatter. 
The  empirical  model  predictions  lie  within  the  scatter  of  the  test  data,  except  for  the  fatigue  test  at 
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loading  level  at  80%  of  ultimate  strength.  All  tests  and  simulations  correspond  to  AS4/3501-6  graphite 
epoxy  composite  from  Shokrieh  and  Lessard  [6]  and  listed  in  Tables  2  and  4. 


Table  4.  Static  material  properties. 


Material  Properties  Strength  Parameters 


Ei 

147  GPa 

X7 

2000  MPa 

e2 

9  GPa 

Xc 

1200  MPa 

e3 

9  GPa 

Yt 

53  MPa 

viz 

0.278 

yc 

204  MPa 

vi  3 

0.278 

ZT 

53  MPa 

^23 

0.40 

Zc 

204  MPa 

G12 

5  GPa 

Sxy 

137  MPa 

Cl3 

5  GPa 

Sxz 

137  MPa 

C23 

5  GPa 

Syz 

42  MPa 

1000  10000  100000  1000000  10000000 
Number  of  Cycles 


10  100  1000  10000  100000  1000000  10000000 
Number  of  Cycles 


Figure  7.  Comparison  of  fatigue  life  of  present  numerical  modeling  and  available  experimental  data 
for  AS4/3501  laminate  subjected  to  tension-tension  fatigue  loading  at  the  frequency  of  10  Hz  and  0.1 
load  ratio,  (a)  [30°]i6  laminate  from  Ref.  [6];  (b)  [90°]8  laminate  from  Ref.  [23];  (c)  [0°/90°/90°]s  from 
Ref.  [23]  and  (d)  [450/-45°]2s  from  Ref.  [33]. 


7.  Open-Hole  Tension  Specimen 

Han  et  al.  [34]  published  experimental  data  for  open-hole  specimen  fatigue  tests  performed  for 
AS4  3501-6  composite  material  laminates.  They  tested  seventy  open-hole-laminate  coupon  specimens 
under  constant  amplitude  tension  and  compression  fatigue  loading.  The  laminate  tested  was  a  24-ply 
quasi-isotropic  [0/  =b  45/90]3s  laminate.  The  tension-tension  fatigue  tests  at  R  =  0.1  were  performed 
for  maximum  stress  levels  of  40%,  50%,  60%,  70%  and  80%  of  the  Notched  Tensile  Strength  (NTS)  of 
the  laminate.  NTS  is  the  ultimate  stress  at  failure  of  the  laminate  under  quasi-static  tension  loading. 
After  one  million  fatigue  cycles  at  40%,  50%,  60%,  70%  and  80%  of  the  NTS,  the  specimens  were  loaded 
to  static  failure  to  obtain  the  Residual  Tensile  Strength  (RTS). 
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The  instantaneous  degradation  scheme  is  known  to  be  mesh  sensitive.  In  fatigue  only  a  small  region 
around  the  hole  sustains  damage  due  to  the  stress  concentration.  However,  to  perform  progressive  failure 
to  get  NTS  or  RTS,  where  damage  leads  to  two  part  failure  of  the  specimen,  the  entire  width  of  the  specimen 
at  the  hole  was  finely  meshed.  The  region  away  from  the  damage  region  where  the  mesh  was  not  refined 
was  modeled  as  elastic  and  did  not  include  the  UMAT  provided  damage  degradation.  This  resulted  in 
a  mesh  with  total  13,600  elements  in  each  ply  and  11,600  elements  in  the  refined  mesh  region  (Figure  8). 

Figure  9  shows  the  stress  strain  response  of  the  laminate  under  quasi-static  loading  in  pristine 
state  and  post  fatigue  tested  state.  The  inset  shows  the  ungripped  portion  of  the  specimen  model  for 
FE  simulations.  Table  5  summarizes  the  comparisons  of  the  predicted  and  experimental  values  of  NTS 
and  RTS  of  the  notched  laminate  specimen. 

The  experimentally  measured  NTS  of  the  specimen  [34]  was  413  MPa.  Progressive  failure  analysis 
prediction  of  NTS  using  the  Hashin's  failure  criteria  and  instantaneous  stiffness  degradation  describer 
earlier  was  405  MPa.  It  is  to  be  noted  the  material  properties  for  AS4  3501-6  used  for  analyses  were  from 
Shokrieh  [6]  as  this  was  not  provided  in  the  reports  of  Han  et  al.  [34].  No  material  property  calibrations 
were  done  to  match  test  results.  Progressive  failure  analysis  of  the  pristine  specimen  shows  that  the 
specimen  experiences  matrix  damage  in  the  90°  and  ±45°  plies  at  40%  of  the  NTS.  This  means  for  fatigue 
loading  above  this  level  there  is  damage  from  the  very  first  cycle.  This  damage  is  higher  when  the  specimen 
is  loaded  at  50  to  80%  of  the  NTS.  Fatigue  damage  for  these  cases  represents  growth  of  the  initial  damage. 


Figure  8.  Mesh  used  for  the  notched  tensile  test  specimen  analyses. 


Figure  9.  Stress-strain  response  of  the  quasi-isotopic  OHT  specimen  in  pristine  state  and  post-fatigue 
tested  state  after  106  fatigue  load  cycles  ( R  =  0.1)  at  different  specified  maximum  stress  levels. 


11 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


J.  Compos.  Sci.  2018,  2,  9 


12  of  21 


Table  5.  Comparison  of  the  residual  strength  of  the  Open  Hole  Tension  (OHT)  specimen  with  test  data  [34]. 


Fatigue  Loading  Maximum  Stress  Level  for 
R  =  0.1,  &  106  cycles 


RTS  Predicted 
(MPa) 


RTS  from  Test  (3  Replicates) 
(MPa) 


40%  of  NTS 
50%  of  NTS 
60%  of  NTS 
70%  of  NTS 
80%  of  NTS 


416 

413 

418 

410 

360 


441, 461, 490 
448, 468,  500 
425, 445, 480 
440, 467, 487 
393, 412, 413 


The  prediction  the  fatigue  residual  strength  RTS  notched  specimens  subjected  to  constant 
amplitude  cyclic  loading  at  40%,  50%,  60%  and  70%  for  106  cycles,  showed  an  increase  in  the  residual 
strength.  The  local  damage  around  the  hole  leads  to  stress  distributions  that  minimize  the  stress 
concentration  at  the  hole.  The  slope  of  the  axial  stress-strain  curve  (Figure  9)  from  progressive  failure 
analyses  of  the  specimens  in  post-fatigued  state  with  simulated  fatigue  damage  shows  that  the  damage 
leads  to  a  relatively  small  stiffness  loss  for  cases  when  the  fatigue  loading  maximum  stress  levels  were 
40%  to  the  70%  of  NTS.  This  is  due  to  the  fatigue  damage  being  primarily  due  to  matrix  failure  and 
damages  being  localized  to  the  region  around  the  hole.  The  stiffness  loss  due  to  fatigue  damage  was 
12%  for  the  specimen  loaded  at  80%  of  the  NTS.  None  of  the  specimens  failed  after  one  million  fatigue 
cycles  at  stress  levels  40-80%  of  NTS.  Post  inspection  of  the  tested  specimens  tested  at  80%  of  the  NTS 
showed  significant  delamination  and  edge  peeling  of  the  plies  [34].  The  delamination  fatigue  damage 
mode  is  not  included  in  our  fatigue  analysis  model. 

During  a  fatigue  loading  simulation,  as  the  strength  degrades  the  failure  index  increases  until  it 
reaches  the  critical  value  of  one,  at  which  point  the  ply  is  considered  to  be  fully  failed.  Figures  10-13 
show  contour  plots  of  the  failure  indices  in  each  ply  for  all  four  failure  modes  after  one  million  fatigue 
cycles  at  40-80%  of  the  NTS.  Failed  elements  are  indicated  by  failure  index  of  1.0.  The  contour  plots  of 
the  failure  indices  show  the  extent  of  failure  and  the  margin  to  failure  at  the  maximum  stress  level 
of  the  fatigue  load  cycle.  Examination  of  these  failure  indices  indicate  that  there  is  significant  fatigue 
damage  in  matrix  tension  mode  in  45°  and  90°  plies  and  matrix  compression  mode  in  the  0°  plies  and 
some  fatigue  damage  due  to  fiber  tension  mode  in  the  0°  plies  at  the  edge  of  the  hole. 

The  fatigue  model  used  here  is  based  on  stiffness  and  strength  reduction  directly  applied  to  the 
engineering  stiffness  constants  and  strengths  that  are  ply  properties.  To  quantify  and  visualize  the 
level  of  damage,  a  measure  of  the  relative  reduction  in  the  stiffness/ strength  parameter  due  to  damage 
Dp,  is  calculated  using  Equation  (20). 
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Figure  10.  Fiber  tension  failure  index  after  106  cycles  in  the  quasi-isotropic  Open  Hole  Tension 
(OHT)  specimen. 
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Figure  11.  Fiber  compression  failure  index  after  106  cycles  in  the  quasi-isotropic  OHT  specimen. 
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Figure  12.  Matrix  tension  failure  index  after  106  cycles  in  the  quasi-isotropic  OHT  specimen. 
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Figure  13.  Matrix  compression  failure  index  after  106  cycles  in  the  quasi-isotropic  OHT  specimen. 
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Figures  14-18  show  the  relative  reduction  of  stiffness  and  strength  properties  due  to  fatigue 
damage  that  exceeded  10%.  The  significant  reductions  (80-100%)  are  observed  in  the  transverse 
to  fiber  (matrix)  direction  moduli  E 2  and  tension  strength  Yj  and  the  plane  shear  moduli  G\2  in 
the  90°  ply  along  the  symmetry  plane  perpendicular  to  the  loading  direction.  A  small  (less  than 
20%)  reduction  in  in-plane  shear  strength  is  observed  in  all  plies.  Although  the  laminate  is  loaded 
in  its  plane,  reduction  of  in-plane  properties  leads  to  transverse  shear  stresses  (captured  by  the  3D 
modeling)  that  induce  damage  and  reduction  in  transverse  shear  strengths  S23.  Although  none  of  the 
specimen  experienced  complete  failure  under  tension-tension  loading  after  a  million  cycles,  radiograph 
images  show  significant  damage.  The  Figures  14-18  show  the  X-ray  radiograph  images  of  the  damage 
accumulated  in  the  specimen  from  Han  et  al.  [34]  at  load  level  of  40-80%  of  the  NTS.  Significant 
damage  is  visible  at  50%,  60%  and  70%  of  NTS,  even  though  the  change  in  residual  strength  is  not 
significantly  high.  At  80%  of  the  NTS  the  damage  state  is  significantly  visible  in  the  specimen  as 
several  modes  of  failure  can  cause  the  damage.  In  addition  to  the  matrix  damage  the  specimen  also 
experienced  delamination  mode  of  failure. 

Figures  19  and  20  show  the  evolution  of  damage  quantified  by  the  relative  reductions  in  the 
transverse  tensile  strength  and  in-plane  shear  stiffness  of  the  ply.  The  progression  shows  that  the 
damage  starts  early  (103  cycles)  and  continues  to  grow  in  a  stable  fashion.  It  is  to  be  noted  that  while 
stiffness  degradations  redistribute  loads  to  other  plies  leading  to  small  decreases  in  stress.  This  has 
the  potential  to  arrest  the  progression  damage.  The  degradation  in  matrix  properties  do  not  seem  to 
play  a  significant  role  on  the  overall  stiffness  or  strength  of  the  notched  specimen.  This  is  because 
the  transverse  (matrix)  stiffness  is  one  order  of  magnitude  lower  than  fiber  direction  stiffness  in  our 
laminate.  However,  quantifying  matrix  damage  is  useful  as  matrix  damage  is  undesirable  as  it  can  lead 
to  delamination.  The  matrix  damage  also  makes  the  composite  prone  to  moisture  intrusion  or  leakage. 
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Figure  14.  Reduction  in  stiffness  E^  in  each  ply  of  the  OHT  specimen  after  106  cycles. 
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Figure  15.  Reduction  in  stiffness  G\2  in  each  ply  of  the  OHT  specimen  after  106  cycles. 


Figure  16.  Reduction  in  strength  Yj  in  each  ply  of  the  OHT  specimen  after  106  cycles. 
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Figure  17.  Reduction  in  strength  S\2  in  each  ply  of  the  OHT  specimen  after  106  cycles. 
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Figure  18.  Reduction  in  strength  S 23  in  each  ply  of  the  OHT  specimen  after  106  cycles. 
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Figure  19.  Reduction  in  strength  Yj  in  90°  plies  of  the  OHT  specimen  for  increasing  fatigue  cycles. 
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Figure  20.  Reduction  in  stiffness  G\2  in  90°  plies  of  the  OHT  specimen  for  increasing  fatigue  cycles. 

The  effect  of  the  interlaminar  failure  and  damage  patterns  are  shown  by  researchers  for  open 
hole  composite  specimens  under  tensile  loading.  Archard  et  al.  [35]  presented  a  discrete  ply  model 
approach  to  experimentally  quantify  damage  under  static  loading  for  a  composite  plate  with  a  hole 
under  tensile  loading.  Nixon-Pearson  et  al.  [36]  presented  an  extensive  experimental  framework  to 
show  the  sequence  of  damage  development  throughout  the  life  of  open-hole  quasi-isotropic  IM7/8552 
carbon  fiber/ epoxy  laminates  loaded  in  tension-tension  fatigue.  Work  of  Nixon-Pearson  et  al.  [36] 
showed  that  there  is  delamination  between  0/45  and  45/90  layers  during  fatigue  cycles.  For  the  open 
hole  tension  specimen,  we  have  depicted  the  reduction  in  the  strength  and  stiffness  properties  in  each 
ply  and  determined  the  residual  strength  of  the  specimen  after  a  million  fatigue  cycles.  The  ultimate 
strength  was  determined  using  Hashin's  failure  criteria  where  right  after  fiber  failure  the  stiffness  E\ 
was  reduced  to  7%  of  its  initial  value  (Table  3).  Thus,  the  effect  of  delamination  to  predict  the  ultimate 
residual  strength  of  the  specimen  can  be  neglected  mostly  due  the  nature  of  the  intralaminar  failure  of 
the  specimen  for  predicting  residual  strength. 

The  methodology  presented  in  this  paper  was  specific  to  the  model  and  material  parameters 
presented  by  Shokrieh  and  Lessard  [17,18].  The  primary  goal  of  this  paper  was  to  develop  and  present 
a  FE  framework  by  using  a  UMAT  subroutine  to  predict  residual  strength  of  the  open  hole  specimen 
and  show  intralaminar  damage  plots  in  each  ply.  The  number  of  parameters  associated  with  the  model 
can  be  significant  sources  of  uncertainty.  The  model  contains  forty  two  parameters  that  is  required 
to  be  calibrated  from  unidirectional  tests.  Thus,  a  separate  study  for  the  uncertainty  quantification 
is  required  for  the  model.  Using  high  fidelity  model  such  as  a  continuum  damage  mechanics  (CDM) 
based  model  has  shown  that  fewer  number  of  such  fatigue  parameters  [37].  The  methodology  (FE 
framework)  can  be  extrapolated  to  any  sets  of  material  provided  the  material  parameters  (Table  2)  for 
static  and  fatigue  tests  are  available  from  experiment. 

Shokrieh's  [6,16,17]  model  develops  the  empirical  S-N  curves,  stiffness  and  strength  degradation 
using  uniaxial  tests  along  each  material  direction  and  uses  these  models  for  multiaxial  fatigue. 
However,  the  model  has  some  limitations.  The  model  does  not  capture  the  correlation  or  interaction 
between  the  degradation  in  stiffness  and  the  residual  strength  parameters  along  different  material 
directions  of  the  ply.  Matrix  cracking  failure  should  cause  degradation  of  both  the  transverse  normal 
and  shear  properties.  In  this  model,  the  transverse  to  fiber  property  degradation  is  only  due  to 
transverse  stress  and  the  in-plane  shear  property  degradation  is  only  due  to  shear  loads.  The  matrix 
cracking  failures  in  the  plies  is  induced  by  both  transverse  to  fiber  normal  stresses  and  shear  stresses. 
The  model  can  be  improved  fitting  the  empirical  models  to  the  evolution  of  damage  variables  with 
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these  damage  variables  used  to  compute  the  constitutive  stiffness  and  strength  properties  following 
a  continuum  damage  approach.  The  models  used  for  fatigue  analysis  of  composites  can  be  categorized 
by  the  fidelity  of  the  model.  In  most  of  the  cases  empirical  models  are  considered  as  lower  fidelity  model 
with  less  computational  time  but  often  having  large  number  of  fatigue  parameters.  Higher  fidelity 
model  includes  the  Continuum  Damage  Mechanics  based  model  and  micro-mechanics  based  models 
that  are  computationally  expensive.  Computationally  cheap  lower  fidelity  models  can  significantly  aid 
in  multi-level  and  multifidelity  modeling  approaches  for  uncertainty  quantification  [38].  Our  revisit 
of  the  Shokrieh  and  Lessard  model  [17,18]  was  for  its  use  as  a  low  fidelity  model  in  a  multifidelity 
framework.  However,  the  inherent  typographical  errors  in  the  original  publication  and  ensuing 
citations  to  the  paper  led  to  large  discrepancies.  In  this  paper,  we  present  the  corrected  version  for 
these  fit  parameters  reported  by  Shokrieh  and  Lessard  [17,18].  As  a  result  of  the  corrections  the  low 
fidelity  model  provides  an  improved  prediction  accuracy  that  can  be  further  refined  and  calibrated 
with  other  high  fidelity  models. 

8.  Conclusions 

This  paper  presents  the  empirical  stiffness  and  strength  degradation  model  of  Shokrieh's  [6], 
its  implementation  as  a  user  material  subroutine  UMAT  in  ABAQUS  FE  software  and  its  validation 
using  test  data  from  literature.  The  results  indicate  this  empirical  model  can  provide  a  good  engineering 
estimate  of  fatigue  life  and  residual  strength  for  preliminary  design  analysis.  The  modification  in 
the  definition  of  the  unified  load  parameter  for  the  transverse  to  the  fiber  direction  improved  the 
prediction  accuracy. 

We  found  that  residual  strength  for  open-hole  tension  specimen  is  mostly  related  to  fatigue 
failure  in  the  fiber  failure  in  0  degree  plies.  Our  developed  user  defined  material  subroutine  (UMAT) 
enabled  us  to  predict  residual  strength  of  the  open  hole  tension  specimen  without  the  inclusion  of 
any  delamination  model.  The  translaminar  failure  is  captured  by  the  degradation  of  the  transverse 
normal  and  shear  moduli  of  the  lamina.  The  inclusion  of  the  cumulative  damage  model  adds  to  the 
capabilities  to  account  for  the  stress  level  changes  due  to  the  failure  during  fatigue  loading. 

FE  based  progressive  failure  analysis  with  instantaneous  degradation  is  highly  mesh  size  sensitive 
and  requires  mesh  convergence  studies.  In  uncertainty  quantification,  the  discretization  errors  must 
therefore  be  included  as  a  source  of  model  uncertainty.  Further,  the  value  of  the  residual  stiffness 
used  after  ply  failure  has  a  significant  influence  on  the  peak  load.  For  example,  changing  the  residual 
stiffness  from  the  values  listed  in  Table  1  to  a  uniform  4%  value  reduces  the  NTS  and  RTS  values  by 
25%.  The  lack  of  correlation  between  degradations  to  transverse  to  fiber  direction  normal  modulus  and 
in  plane  shear  due  to  matrix  damage  is  also  a  source  of  model  error.  These  two  are  the  most  significant 
source  of  epistemic  uncertainty  that  needs  higher  fidelity  models  and  test  data  to  correct. 
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